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Of the various mechanisms put forward in recent years to explain long range tropospheric 
radio propagation the most important seem to be scattering from atmospheric turbulence 
and partial reflexion from high level subsidence inversion layers. In this paper the writer 
extends earlier theory appropriate to the case of a single elevated inversion layer to cover the 
case of multiple layers. In some cases it is found that the theory and formulas can be applied, 
with very little modification, to the latter. The special cases in which this can be done are 
worked out in detail and it is found that, if certain conditions are satisfied, several weak high 
level inversions can produce a similar effect on the propagation to a single strong inversion. 



1. Introduction 

Of the various mechanisms put forward in recent years to explain long range tropospheric 
propagation the most important seem to be scattering from atmospheric turbulence and partial 
reflect ion from high level subsidence inversions. In 1952 the writer, after first reviewing earlier 
experimental work, showed theoretically that a well developed inversion could cause the phe- 
nomenon. In a three part paper in 1955, after first reviewing the major theories which had 
been put forward in the interim, he considered in greater detail the special characteristics of 
the field under a high level subsidence inversion (part II). It is now fairly generally agreed 
that the major mechanisms responsible for the long distance fields in tropospheric propagation 
are scattering from atmospheric turbulence and partial reflection from high level subsidence 
inversions. 

Refractometer measurements have shown that several minor high level inversions can 
coexist at different levels and that this situation is probably present for most of the time. It 
would therefore seem desirable to make the previous investigations complete by extending the 
theory to take account of the case when several such inversions are present. This is what will 
be attempted in the present paper. 

2. Propagation Under Several High Level Inversions 

It is possible to show [Northover, 1952] that, provided the distance of the field point from 
the exciting source is small compared with the radius of the earth, the Hertzian function n 
representing electromagnetic disturbances excited by a horizontal electric dipole (or system of 
dipoles) in a stratified atmosphere outside a spherical earth is expressible in the form 

-Ti^n(kr)P n (cose) (1) 

r in) 

where ^l n (kr) is a solution of the equation 



1 Contribution from Department of Mathematics, Carleton University, Ottawa, Canada. 
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and that this representation is valid whether the dielectric constant k or the refractive index jn 
is a single analytic function of height (as for a surface inversion) or has different forms for 
different regions (as when high level inversions are present). Further, it appears that, subject 
to the same restriction on range, the condition at the surface of a perfectly conducting earth 
can be taken as II— 0. 

Consider the case when several high level inversions are present. Let there be p of them 
and let the boundaries of the jth be given by $,=y h t=0jfa>yi); let also 



^o(=x)<y 1 <z l <y 2 <^z 2 



<Vp<z p < • 



(3) 



Let the dielectric constant below the lowest 2 be k 2 and above the highest be unity; let also 
k=k 2 in 2i<£<2/2, K=K3in ^<^<y 3 . ■ . k=kj in ^-i<£<^ and k=k p in z,-i<£<iyp. 

These values of k refer to the space between the layers; throughout a given layer k is 
supposed to change smoothly, e.g., through the second layer k changes smoothly from k 2 at 
the bottom (%=y 2 ) to /c 3 at the top (£=z 2 ). See figure. 




e= y , 



Remembering (1) we assume that the secondary field below all inversions (i.e., in the 
space £<£<2/i) is of the form 



where 



-^S(2»+l){2^f,ttV^)+e r .i?.ttV^)}i > .(ooB^ 

(1 \ 1/2 /l Y /2 

-TXj H%\ V2 {X)) Vn(x)=^7rX) #4+1/20*0 

and we are assuming the time factor e ioit , and the primary contribution from the dipole is 



(4) 
(4A) 



and 



where 



-^S(2w+l)f„(6i)^(fV^)P«(cos 6), when £V*i<&i 



-^S(2n+l)^(6 1 )f B (f V*i)P»(cos 0), when f V^>6i 



rbi 



b 1 =2T(R+H)^i Kl /\;2Mx) = Ux) + Vn(x). 



(5) 



(5A) 



2 The normal slow lapse rate of y. will be approximately allowed for by taking an "effective" value for the earth's radius. This is a reasonable 
procedure if the inversions are not too high [cf Norton, 1959; Bremmer, 1960; Wait, 1959]. 
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The constant f3 is proportional to the square root of the power radiated: (see reference 1, 
page 111), J? and i?are, respectively, the earth's radius and antenna height. 
Let the total field in the space s./-i<Cf<J/./> i-e., the region where k=kj be 

-^Z!(2^+l){^ ) f w (^)+5f ifc(*Vfd }P n (cos 9)... (6) 

and within the jth inversion layer (y ; <Cf <C^) let it be 

-^S(2n + l){C« , <a)+Z?« ) ^'(|)} J P re (cos 9) (7) 

where u^(^), v^ (£) are any pair of fundamental solutions of (2) within the jth layer {yj<i 
The field above the complete set of inversion layers is of the form 

-^E(2n+l)E n U0Pn(.cos 6)- (8) 

It will be useful to note here that, in view of the continuity of tin 4 rising and descending field 
contributions at the lower boundary of the first layer, we have the formulas 

B™=G n (9) 

for &i<C£<Cyi) am l ll l s0 ) that 

A<>+»=E H ; B*+»=0 . . . (10) 

in view of the radiation condition at infinity. 

3. Boundary Conditions 

These conditions, expressing the continuity of the field and of its normal derivative at the 
interfaces of the successive inversion layers, read as follows: 

At£=2,_, 

where j=2, 3, . . ., p. 

C^ 1, tf- 1 1y^l)+5^ 1 ^-^!/^,)=4^ 1) f,(!/^l^- 1 )+^^^»(!/^ 1 ^- 1 ) ) , 

+^-"^-^0} 

wherej=3,4, . . .,p+l. 
If we let 

B<P/A!P=-H t . . . (13) 



(11) 



we find, after somewhat involved algebra 






(14) 
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where 

(?,_ 1 =iv«- i) rJ^_ 1 VK^)rn(^-iV^)+i - 1> f;(^- 1 V^)f;(^-iV^) 

-M 1 °- i ^;(^_iV^i)r»(sy-iV^)--M 2 w ~ 1) '7 n (^-iV^)f;(2y-iV^) (i4a> 

the bar means that we interchange the f and 77 functionality symbols in these formulas, and 

Afl»={ttf>(«,) } Wtod-uPM {vkHzj) } ' (14b) 

Mf=J^ {«£»(«/) {t>« (2/,) }'-{«# (y,) } 't>« ( 2 ,) } 

^" -^ [{«-"' (*i) }' K^W-W' (</,)}' {^ 0) (^)}'J. (14B) 

In cases where the earth behaves effectively as a perfect conductor the surface condition is, 
as already mentioned, II =0, by (4) and (5) this amounts to 

JK+lUbi) \Ux^)+^G n +^Ubi)\vn(^0=O. (15) 

Hence, by some rather involved algebra, it is possible to derive an expression for the total 
Hertzian function at points below the lowest inversion similar in form to that developed in 
the writer's 1952 paper, namely, 

n=-2^2 (2 ^ 1) Ux!)-h!v!S) {f-^'-^-^^f.fe)}^^*) (16) 

in which 

£i=fV^i5 Xi=x4k\\ x=2tR/\ (16a) 

and Hi is the value of Hj when^=l : it is expressible in terms of F n and G n by (9). 

We have now arrived at the mathematical solution of the problem, for, since H x can now 
be determined (see the argument immediately following), then (9) and (15) now amount to 
two simultaneous linear equations for F n and G n . Determination of H x : 
We have, from (10) 

■*-* n 

A (p+1) 

=0by(10) (17) 

we have 

Hp=Qp/P P (17a) 

and then all the H's from H v „ x down to H x are given in succession by (14). 
4. Approximate Evaluation of the Field 

The field series (16) cannot, as it stands, be used to discuss the behavior of the field 
below the inversions because of its slowness of convergence. To make further progress, it 
is necessary to transform it into a more rapidly convergent series (series of residues) by complex 
^contour integration (Watson's transformation), as was done in the writer's 1952 paper on the 



TT 

tJ-p + 1 A (p+1) 

J - ± n 



single high level inversion. As in that paper, the contribution of the imaginary axis to the 
complex contour integral involved is zero, and it appears that the series of residues giving 
the field is of the same general form as that obtained therein. It is as follows: 

7T^ ,« fp,- -1/2(6)^-1/2(31) — £,,-1/2(3:1)^-1/2(6) 

l ~rb U V} ~ 2i-{ Vv .^ /2 (x l )} 2 [dH l /ds] {s=Vj) 

X{r^-i/2(»i)^-i/2(|i)— 17^-1/2(^1) f^-i /2 (?i)}sec vpP 9r ui(— cos 6) (18) 

where the summation is over the roots vj of the equation 

{ s -m(xi)/ris-i/2(xi) = H 1 (19) 

considered as an equation in s. 

5. Thin Layers — Meter Wave Propagation 

It can be shown that, as in the case of a single inversion layer, the high level inversions 
have their maximum effect upon the propagation when they are sufficiently thin to act effec- 
tively as dielectric discontinuities. The condition that this may be so is that I^2p must be 
small compared with unity for each inversion, where I=2tt/\, p=--h/R being the layer thick- 
ness and h the inversion height. We shall, therefore, confine ourselves to a discussion of this 
case only. 

When the layer thickness may be neglected, the preceding results simplify; we have 

M[*-» ^-" = -4^7,:^ (20) 

and since only the mutual ratios of L u ~ l \ Mi u_1) , M 2 {j ~ l) and N u ~ l) are significant, we may 
take 

Qj-i = Vfy-i K-1/2 (jO-iVty-0 ?*-i/2 (%-iV^) — V^^-1/2 \yj-i-\Kj-i) fJ-1/2 (zj-i-Jty) 

Pj-1= Vty-1 ^-1/2 (Vj-l^Kj-l) fp-1/2 (^-lV^)""V^^-l/2 (2/y-lV^-l) ri-1/2 (%-lV^) 
^_l=V^-l^-l/2 GA-iV^-i) ^-1/2 (zj-I-Jk}) — -V^^-l/2 (2/y-lV^-l) *7*-l/2 (zj-itJkj) 
Pj-1= V^l fi-1/2 (j/i-1 Vty-0 ^7,-1/2 (%-l Vty) """" Vfy J"* -1/2 (2//-1 V^^-l) ^C — 1/2 (Sj-1 Vty) (21) 

It is possible to make further progress only when the heights of the high level inversions 
are large compared with R 1/3 \ 2/3 . The inversions are then to be supposed as situated, above a 
well-known critical height which is of great importance in all diffraction theories for propa- 
gation along a sphere; asymptotic approximations may be applied there to all height gain 
functions. Fortunately, this condition is usually satisfied in practice, even for meter length 
VHF waves. It is then possible to show that the only roots v which could possibly give im- 
portant contributions to the field series (18) by reason of smallness of — Im (v/x 1/z ) fall into the 
following types: 

(A) those near 3 £ =x 

(B) those near £=?/,■ 0=1, 2, . . . , p) 

(A) corresponds to waves nearly grazing the horizontal level at the receiver and (B) to waves 
nearly grazing the lower boundaries of the successive inversion layers. 

It can further be shown that the contribution of roots of the group (B) to the field series 
is unimportant compared with that of group (A). Although limitations of space forbid the 

3 By "$ near z" we mean ||-z| - Q(xM). 
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inclusion here of the exact analysis, it may be remarked that the reason for this is to be found, 
not in any wide discrepancies in the magnitudes of their imaginary parts, but in the fact that 

{4 / b n \ 3/ 

-■yj2x (e—yj V 

i.e., of order Exp (-^2xp) ,2 \ 
where 

Pj=(yj—x)/x, by virtue of the factor {r} v _i l2 (x\)} 2 {we have written e=(v—x)/x in the above}. 

This makes the leading terms of the series (18) which correspond to values of v in this 
group exponentially small. 

We shall therefore confine ourselves to consideration of values of v for which \v— x x \ is 
0(# 1/3 ). 4 For such values, we may apply well known asymptotic properties of Bessel functions 
(cf Watson's treatise or the writer's previous paper, 1) and we obtain, from (21) 5 

+d '" l}1/2 ] * {2(pJ-<)V» Exp {-iVftKft-.-.)^} 

K7^P,- 1 M[{2(^- 1 -e) + d,}^+( K ^) 1/2 {2(^- 1 -e) + d y _ 1 } 1/2 ]- {2(y; _ i 1 _ f)} i/2 



+ d. 



«7 1/2 Py-i=-*[{2(^- 1 -e) + d ; }^+('^i) I/2 {2(^_ 1 -e) + d,_ 1 }^]. 



{2(^-x-0} 1/2 



(22) 



where 



kj=l + dj, etc. (22A) 

Using these approximations in (14), we obtain, since vjx may be taken as effectively unity 



Hi-i= 






}-t— ~ 



[{2( 2> ,_ 1 - e )+d,} 1 / 2 +{2(^_ 1 - e )+d,_ 1 } 1 « ! ]-iiy i [{2(^_ 1 -6)+d,}^ 

-{2(^_ 1 - 4 ) + V 1 } 1 / 2 ]Exp/|v2a ; i( 2 ) y -6)^\ 



] 



(23) 



It transpires that Im (e) {which is always negative, so that Pj-i— cis situated in the first quad- 
rant} is never very small compared with x~ 2/3 , but pP>x~ 2/d and hence 



{- 



Exp ^-|V2^(^-!-e) 3/2 H>1; 



»}l 



Exp ^xi{ Vi ^- ( y- 



•CI. 



Again, it appears that H, is never large compared with unity. The denominator may then be 
effectively replaced by 

{20p,_ 1 -«)+d,} 1 'H-{2fo_,-«)+d,_i}»'» 



4 See also note 2. 

5 Since zpf 2 »l- See note 3. 
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and then (22) and (23) give 

Hj.^kj^+Hj (24) 

which yields at once 

H 1 =k 1 +k 2 + ...+k p (25) 

where 

hQj/Pj, by (22). (25a) 

But kj is the value which H x would have in the field series (18) if the jth. inversion only were 
present. We obtain, therefore, the following important result: 

Addition Theorem. "When the inversion layers are at heights large, or moderately large 
compared with i? 1/3 X 2/3 , and when also they are sufficiently thin to act effectively as dielectric 
discontinuities, then the value of Hi to be used in the field series (18) and in the mode equation 

(19), is 2i2kj where kj is the value which Hi would have in (18) if thejth inversion only were 

present." 

6. Discussion of Results 

Even when the conditions of the above theorem are fulfilled, the resulting equation (19) 
for v is, unfortunately, still too complicated to allow of much further progress being made. 

There is, however, a special case of the theorem in which the characteristic properties of 
the propagation can be clearly perceived. It occurs when the layers are relatively close 
together. 6 As far as is known, this case is likely to occur in practice. It might, therefore, be 
useful to conclude by considering it. 

7. Case When Inversions Are Close Together 

In practice the dielectric constant differences 7 bj are all small compared with pj. Since 
by hypothesis xpf'^l and e is to be 0(x~ 2/3 ), (25A) gives the following approximate formula 
for kj: 



*> = -*I^M~i^ (2, '~ €)3/2 } 



(26) 



When the inversions are close together, we may take pj in (26) to be effectively independent 
of j; that is, we may take pj equal to an "effective" value p E . It would be appropriate to call 
Rp E the "cluster height" or "group height" of the system and this height would not differ much 
from the heights of any of the inversions in the case we now have in mind. Thus we have, 
in this case 



H ^-^) E ^{-h r2xi{ ^- e)3/2 } 



(27) 



where A is the total drop in dielectric constant across all the inversion layers. 

Hence, in this case, the inversions behave like a single inversion layer across which the 
refractive index drop is equal to the algebraic sum of all the differences of refractive index 
across the separate inversions of the group. 



6 I.e., their widest separation is small compared with the height of the lowest. 

7 These are twice the corresponding refractive index differences. 
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8. Conclusion 

In the writer's 1952 paper, it was found that a sharp high level inversion whose dielectric 
constant difference was about 2.58 X10~ 5 could cause considerable anomalous propagation of 
meter radio waves. Inversions of this size do not, perhaps, often occur, but under settled weather 
conditions it is more likely that sharp high level inversions of much smaller dielectric constant 
difference would occur much more frequently. If several such small inversions occurred simul- 
taneously {e.g., four such that dj— d J+ i = 0.6X10" 5 would do}, the above theory shows that 
the propagation would be effectively the same as that characteristic of a single large inver- 
sion. A mechanism has therefore been found which is capable of explaining the long distance 
fields without having to invoke the existence of a single strong elevated high level inversion. 

Note 1. The influence of scattering from atmospheric turbulence. The present paper has shown that 
when several high level inversions are simultaneously present their effect on propagation can be equivalent 
to that of a single well-marked high level inversion aloft, and Saxton [1951] in his experimental work on 9 and 
45 Mc/s waves has found that the effect of scattering is, in general, likely to be unimportant compared with reflec- 
tion effects from (the equivalent of) a single inversion aloft, at any rate over the path considered by him. On 
the other hand, experiments elsewhere seem to indicate that scattering from turbulence can be at least as 
important a factor in the propagation [Northover, 1955]. 

The basic meteorological question is whether departures from the standard atmosphere structure, which 
are significant from the radio propagation point of view, occur in a stratified form (which may be viewed as 
anisotropic turbulence), or as isotropic turbulence. The present paper has attempted to make the theory 
sufficiently general to cover the most important possibilities of the first case. The best type of radio meas- 
urement for deciding between these alternatives appears to be space diversity measurement in mutually per- 
pendicular dimensions. The difference in the angular scattering characteristics of the two models could be 
ascertained, in theory, from such measurements. Experimental work along these lines has been carried out 
and results which have been obtained so far seem to favor the stratified model. 

Note 2. On the location of the roots v. When several high level inversions are present simultaneously 
we have found that there is one case only which is tractable mathematically, namely, that in which the inver- 
sions behave like a single high level inversion. This, therefore, is the only case which has been worked out 
in detail in the text and the disposition of v for it is known from the writer's 1952 paper dealing with a single 
inversion layer. It was shown therein that the significant values of v were those "near" x\ 

{i.e., 1^^1 = 0(^/3)}. 

Note 3. Asymptotic approximations to Hankel-Nielson functions. These are summarized on page 116 
of the writer's 1952 paper. Using them we obtain, since 

xp^ 2 » l;d,--i<<P/-i;z€ 8 ' 2 = (1) 

ir / r \ ( 1 1-1/4 |2/—/ 1 \ 3/2 ) 

2 1/4 e — $v-m\Vi-\SKi-\)~\ flf-i — €+^d*_i | exp | —~^2ix ( py-i-€-f ^^-i ) \ 
2 1/ %-xr:-i/ 2 (l/,--iV^T)-|pf-i-e+|d y _ 1 1 1/4 V2 exp j -\^2ix (pi-i-e+^d^y^ 
2 1/ ^-x^-i/ 2 (^-iV^)-(p,-i+^-€+|d J )" 1/4 exp { -IV^^^-h^-c+id,) 372 } 

2^-xf;-i/2U-iV^)^^^ 

in which ?-,• is the thickness of thejth inversion. The corresponding expressions in r), V are obtained from the 
above by changing the sign of the coefficients of i in the above. 

This enables us to obtain (22) from (21) neglecting the term tj-i/R and the terms d,_i, d,- in comparison 
with p^ i in the exponentials. 
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